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Abstract — The internal flow field of a square container and a 
circular container was simulated by finite element software. The 
mixing result of fluid in the container and the influence on 
growth of the floes during coagulation mixing were analyzed 
from the following four aspects: average turbulent dissipation 
rate, average speed, average turbulent kinetic energy and the 
scale of micro vortex. The result shows that the distribution of 
velocity field, turbulent kinetic energy field and turbulent 
dissipation energy field is more balanced in square container 
and leads to the formation of more X scale of micro vortex. The 
hydrodynamic conditions are more favorable for the growth of 
the floes. 

Index Terms — coagulation mixing; flocculation; scale of 
micro vortex; numerical simulation 

I. Introduction 

In the process of water treatment, the growth of the floes will 
affect the efficiency of water treatment directly [1]. There are 
many factors such as the shape of the container, mixing time, 
mixing speed, the size and shape of the mixing paddle that 
affect the growth of the floes [2, 3]. In this paper, the 
three-dimensional flow field in the container was calculated 
and simulated by FLUENT software. The effect of circular 
and square container on the growth of the floes was studied 
while the mixing speed, mixing time, and mixing paddle 
shape remain unchanged. The calculated turbulent dissipation 
rate and scale of micro vortex were used as the evaluation 
index of the flocculation effect [4] and provide some 
theoretical guidance for the practical engineering 

II. Numerical simulation research of mixing 

Multi-reference system model ( MRS ) [5] was adopted in 

the numerical simulation. The specific research process is as 
follow: a fluid rotating reference frame which including the 
mixing paddle and the surrounding fluid was defined, and 
then a stationary coordinate system which was outside the 
area of the mixing paddle was defined in the container. The 
shape is shown in Fig 1 and Fig 2 (The dotted line means 
section of the two reference frame). Assuming that the section 
of the two reference frame is steady flow, that means the 
speed (absolute speed) of the section for each reference frame 
must be the same, and the grid does not move. 



Fig.l Sketch Map of mixing paddle in the circular 

container 



Fig.2 Sketch Map of mixing paddle in the square 

container 

III. Establish mesh 

The hexahedral mesh was chosen in the simulation model 
and the minimum interior angle of each mesh cell was above 
50 DEG. All determinant (2*2*2) options in mesh cell were 
above 0.7. The mesh models are shown in Fig 3. The k-s 
equation model was used in the calculation model, the 
boundary condition was the whole wall function. 



Fig.3 Mesh model of circular container and square 

container 
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Internal flow Numerical Simulation Analysis of Coagulation Mixing 


IV. Mixing simulation comparison between circular and square container 
The mixing speed was set to 65 rpm. 

A. Mixing paddle 1 ( 40><40mm ) 


Simulation cloud image of mixing in the circular container : 
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Contours of Turtxient Dissipation Rate (Epsilon) (m2/s3) 


Contours of Velocity Magnitude irrVs) 


Contours of Turbulent Kinetic Energy (k) (m2/s2) 


a. the turbulent dissipation rate b. the turbulent velocity 


c. the turbulent kinetic energy 


Fig.4 The distribution of three indexes in circular container 


Simulation cloud image of mixing in the square container : 


I 5.61e-02 
5 33»-02 
5 05*02 
4.77*42 
4 49*42 
421*42 
393*42 
3 65*42 
3.37*42 
3 09* 02 
2 81*42 
2.53*42 
2 25*42 
1 .97*42 
1 69*42 
141*42 
m 113*42 
845*43 
5 65*43 
2 85*43 
4 93*45 




1 2.17*41 
2 06*41 
1 96*41 
1 85*41 
— 1 74*41 
■ 1 63*41 

1 52*41 
141*41 
1 30*41 
1 20*41 
1 09*41 
9 78*42 
8 69*42 
7 81*42 
6 52*42 
543*42 
4 35*42 

y — 326*42 
217*42 
1 09*42 
0 00*^00 




I I 27*43 
121*03 
1.14*43 
1.08*43 
1 .02*43 
9.58*44 
8.96*44 
8.34*44 


7.72*44 

7.10*44 

6.48*44 

5.86*44 

524*44 


4.62*44 



4.00*44 

3.38*44 

2.75*44 

2.13*44 

1.51*44 

8.93*45 

2.72*45 



Contours of Turt>utent Dissipation Rate (Epsilon) (m2/s3> 


Contours of Velocity Magnitude (rrVs) 


Contois of Turtxjlenl Kinetic Enerav (K) (m2/s2) 


a. the turbulent dissipation rate b. the turbulent velocity c. the turbulent kinetic energy 

Fig.5 The distribution of three indexes in square container 


B. Mixing paddle 2 ( 50 x 40mm ) 


Simulation cloud image of mixing in the circular container : 




1 40*4)1 
1 33*4)1 
1 26*411 
1 19*4)1 
1 12*01 

1 05*411 

6 80*4)2 
910*4)2 
841e4)2 

7 71*4)2 
781*4)2 
6 32*4)2 
5 62*4)2 
4 93*4)2 

4 23*02 
3 53*412 

2 84* 02 
2 14*02 
1 44*4)2 
747e4)3 

5 04e4M 




1 283 * 4)1 
268*411 
254 * 4)1 
240 * 4)1 
226 * 4)1 
2 . 12 * 4)1 
1 98 * 4)1 
1 84 * 4)1 
1 70 * 4)1 
1 56*01 
141 * 4)1 
127 * 4)1 
1 . 13 * 4)1 
9 89 * 4)2 
848 * 4)2 
7 06 * 4)2 
666 * 4)2 
, _ 424 * 4)2 

I ™ 1 283 * 4)2 
141 * 4)2 
0 00**00 




1 293*4)3 
2.79*4)3 
2 65*4)3 
251*4)3 
2 38*4)3 
224*4)3 
2.10*4)3 
196 * 4)3 
1 82 * 4)3 
1 .68*4)3 
1 . 54 * 4)3 
1 40*4)3 
126*4)3 
1.12*4)3 
9 86*4)4 
8 47*4)4 
7 08*4)4 
. _. 5 69*4)4 
H 430*4)4 
291*4)4 
• 1 52*4)4 



Contours of Turtxient Dissipation Rale (Epsilon) (m2/s3) 


Contours of Velocity Magnitude (m/s) 


Contours of Turbulent Kinetic Energy (k) (m2/s2) 


a. the turbulent dissipation rate b. the turbulent velocity c. the turbulent kinetic energy 

Fig.6 The distribution of three indexes in circular container 


Simulation cloud image of mixing in the square container : 
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a. the turbulent dissipation rate b. the turbulent velocity c. the turbulent kinetic energy 

Fig.7 The distribution of three indexes in square container 

Mixing paddle 3 ( 50 x 50mm ) 


Simulation cloud image of mixing in the circular container 
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Contours of Turtulenl Dissipation Rate (Epsilon) (m2/s3) 
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a. the turbulent dissipation rate b. the turbulent velocity c. the turbulent kinetic energy 

Fig.8 The distribution of three indexes in circular container 

Simulation cloud image of mixing in the square container : 
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a. the turbulent dissipation rate b. the turbulent velocity c. the turbulent kinetic energy 

Fig.9 The distribution of three indexes in square container 

D. Mixing paddle 4 ( 50 x 60mm ) 


Simulation cloud image of mixing in the circular container 
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Fig. 10 The distribution of three indexes in circular container 

Simulation cloud image of mixing in the square container : 
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Fig. 11 The distribution of three indexes in square container 

E. Mixing paddle 5 ( 60 x 40mm ) 


Simulation cloud image of mixing in the circular container 
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a. the turbulent dissipation rate b. the turbulent velocity c. the turbulent kinetic energy 

Fig. 12 The distribution of three indexes in circular container 

Simulation cloud image of mixing in the square container : 
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Contours of Turtxient Dissipation Rate (Epsilon) (m2/s3) 
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a. the turbulent dissipation rate 


Contours of Turbulent Kinetic Energy (k) (m2/s2) 

c. the turbulent kinetic energy 


b. the turbulent velocity 
Fig. 13 The distribution of three indexes in square container 

Notice: the above cloud images are the cross section which is on Y axis of the container and 50mm from the bottom. 


Comparing with the cloud images above, it can be found 
whether mixing in circular or square container, the turbulent 
dissipation rate field, turbulent velocity field and the turbulent 
kinetic energy field are concentrated around the mixing 
paddle, and gradually expand outward, but there is a large 
difference between the two containers in expansion degree. 
Comparing with each set of simulation one by one, the index 
distribution is more balanced in the square container under the 
same mixing conditions, the most obvious difference exist in 
the turbulent dissipation rate. 

For a more intuitive analysis of the various indexes, line 
charts of the turbulent dissipation rate, average speed and the 
turbulent kinetic energy which were solved under the mixing 
conditions by the FLUENT software were drew for 
comparison. The results are shown in the following three 
figures. 


Fig. 14 Comparison of average turbulent dissipation rate 
between circular and square containers 
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Fig. 15 Comparison of average speed between circular 

and square containers 
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Fig. 16 Comparison of average turbulent kinetic energy 
between circular and square containers 
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It can be seen from the chart above, the average turbulent 
dissipation rate, average speed and average turbulent kinetic 
energy in circular containers are larger than the average 
values in square container under the same mixing condition. It 
can be seen obviously that the turbulence intensity in the 
circular container is larger than the turbulence intensity in the 
square container, from the perspective of the floe collision 
theory [6], the floe collision will dissipate more energy in the 
circular container than the square container, and it also 
generate greater shear force which will break floes; from the 
perspective of the scale of micro vortex theory [7], the range 
of the vortex scale is wider in circular container, but when the 
ratio of scale of micro vortex which is similar to the floes gets 
smaller, the growth of floes gets more difficult. 


V. Comparison of the scale of micro vortex 


By the formula /l = {v* / 


1/4 


the data which were 


obtained by numerical simulation are used to draw a line chart 
of the scale of micro vortex in two different containers. As 
shown in Fig 17. 



Fig. 17 The line chart of the scale of micro vortex in 
different containers and by different mixing paddle 


It can be seen from the chart, whether in a circular container 
or a square container, the tendency of the scale of micro 
vortex which changes with the mixing paddle is same as the 
tendency of floe changing (Fig 17). In addition, under the 
same condition, the scale of micro vortex in a circular 
container is smaller than the square container. Thus, the 
distribution of the turbulent dissipation rate is more balanced 
in the square container, the scale of micro vortex in the fluid is 
bigger than the circular container. The conclusion can be 
reached that the effect of coagulation mixing in the square 
container is better than the circular container by the scale of 
micro vortex collision theory [8]. 
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VI. CONCLUSION 

The internal flow field of a square container and a circular 
container were simulated by the finite element software. The 
result shows that the distribution of velocity field, turbulent 
kinetic energy field and turbulent dissipation energy field is 
more balanced in square container under the same condition. 
The hydrodynamic conditions are more favorable to the 
formation of the flocculation. Meanwhile the square container 
generates more X scale micro vortex which will lead to the 
formation of the coagulation and make a contribution to water 
treatment. 
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